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Computational  Approach  To  Characterize  Causative  Factors 
and  Molecular  Indicators  of  Chronic  Wound  Inflammation 


Sridevi  Nagaraja,  Anders  Wallqvist,  Jaques  Reifman,  and  Alexander  Y.  Mitrophanov 

Chronic  inflammation  is  rapidly  becoming  recognized  as  a  key  contributor  to  numerous  pathologies.  Despite  detailed  investigations, 
understanding  of  the  molecular  mechanisms  regulating  inflammation  is  incomplete.  Knowledge  of  such  critical  regulatory  processes 
and  informative  indicators  of  chronic  inflammation  is  necessary  for  efficacious  therapeutic  interventions  and  diagnostic  support  to 
clinicians.  We  used  a  computational  modeling  approach  to  elucidate  the  critical  factors  responsible  for  chronic  inflammation  and  to 
identify  robust  molecular  indicators  of  chronic  inflammatory  conditions.  Our  kinetic  model  successfully  captured  experimentally 
observed  cell  and  cytokine  dynamics  for  both  acute  and  chronic  inflammatory  responses.  Using  sensitivity  analysis,  we  identified 
macrophage  influx  and  efflux  rate  modulation  as  the  strongest  inducing  factor  of  chronic  inflammation  for  a  wide  range  of  scenarios. 
Moreover,  our  model  predicted  that,  among  all  major  inflammatory  mediators,  IL-6,  TGF-P,  and  PDGF  may  generally  be 
considered  the  most  sensitive  and  robust  indicators  of  chronic  inflammation,  which  is  supported  by  existing,  but  limited,  exper¬ 
imental  evidence.  The  Journal  of  Immunology,  2014,  192:  1824-1834. 


Inflammation  is  an  essential,  nonspecific  innate  immune  re¬ 
sponse  that  facilitates  survival  during  infection,  injury,  and 
disease  (1,2).  In  the  absence  of  a  persistent  initiator,  a  well- 
balanced  inflammatory  response  usually  resolves  in  ~2  wk  (3,  4). 
Inflammation  is  normally  reduced  when  infiltrated  leukocytes  are 
eliminated  from  the  inflamed  site  and  the  tissue  populations  of 
macrophages  and  lymphocytes  return  to  their  normal,  prein¬ 
flammation  numbers  and  phenotypes  (2,  3,  5).  The  inflammatory 
response  is  coordinated  by  a  group  of  cell-derived  molecular 
species  known  as  cytokines,  chemokines,  and  growth  factors. 
Most  cytokines  are  multifunctional  and  are  involved  in  an  ex¬ 
tensive  network  that  carries  out  positive  and  negative  regulation  of 
cell  activation  and  behavior  (6-8).  Dysregulation  of  the  inflam¬ 
matory  response  may  lead  to  prolonged  arrest  in  an  inflamed  state, 
possibly  resulting  in  host  tissue  damage  and  pathological  chronic 
inflammation. 

Chronic  inflammation  can  be  generally  defined  as  a  physiolog¬ 
ical  state  that  is  characterized  by  a  prolonged  and  heightened  in¬ 
flammatory  response  (9,  10).  It  has  been  associated  with  various 
autoimmune,  infectious,  and  neoplastic  pathologies,  including 
tumorigenesis,  CNS  disorders,  diabetes,  psoriasis,  atherosclerosis, 
rheumatoid  arthritis,  asthma,  and  chronic  wounds,  among  many 
others  (2,  10,  11).  Despite  its  central  role  in  disease,  the  mecha- 
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nistic  details  underlying  the  initiation  and  progression  of  chronic 
inflammation  remain  largely  unresolved.  In  addition,  reliable 
molecular  indicators  for  the  early  detection  of  chronic  inflam¬ 
mation  have  not  been  identified  because  of  a  lack  of  relevant  data 
for  human  subjects.  Chronic  inflammatory  pathologies  have  many 
manifestations,  such  as  abnormal  apoptotic  neutrophil  loading  in 
diabetic  ulcers  (12)  and  increased  presence  of  classically  activated 
macrophages  in  chronic  ischemic  wounds  (13).  Yet,  the  most 
conspicuous  and  often-reported  feature  of  many  chronic  (or 
delayed)  inflammatory  scenarios  is  the  persistent  presence  of  el¬ 
evated  levels  of  pro-inflammatory  cytokines,  such  as  TNF-a,  IL- 
1(3,  and  IL-6  (12,  14-17). 

Animal  model  studies  of  chronic  inflammation  in  wounds  and 
diseases,  such  as  obesity  and  atherosclerosis,  have  elucidated 
gene-  and  protein-level  differences  between  normal  and  chronic 
inflammatory  responses  (12,  13).  Yet,  it  has  not  been  firmly 
established  whether  the  expression  of  the  same  genes  is  altered 
in  humans  during  disease,  questioning  the  validity  of  such  in¬ 
dicators.  Any  human  study  of  chronic  inflammatory  diseases 
takes  place  well  after  the  disease  has  already  progressed  and, 
thus,  fails  to  provide  insights  into  the  factors  triggering  the 
chronic  state  (9,  18-22).  Although  chronic  inflammation  can  be 
present  in  many  pathological  conditions,  it  is  not  causative  to 
disease  unless  exacerbated  by  additional  genetic  and  environ¬ 
mental  factors  (11).  Therefore,  attempts  to  alleviate  chronic 
inflammation  are  confined  to  cases  in  which  a  disease  has  re¬ 
sulted,  at  which  point  the  treatment  is  focused  on  addressing 
the  symptoms  of  inflammation  as  opposed  to  its  causative  fac¬ 
tors.  As  a  result,  many  basic,  mechanistic  questions  regarding 
the  chronic  state  of  the  inflammatory  process  have  remained 
unanswered. 

Because  of  the  inherent  complexity  of  the  inflammatory  process 
and  the  aforementioned  limitations  of  conventional  approaches, 
these  questions  are  difficult  to  answer  via  traditional  experimen¬ 
tation  alone.  A  complementary  in  silico  approach  has  the  unique 
advantage  of  providing  focused  and  time-efficient  integration  and 
analysis  of  the  available  literature  data,  with  the  capability  of 
generating  experimentally  testable  hypotheses  to  expedite  the  in¬ 
vestigative  process.  Although  a  number  of  mathematical  models 
recently  have  been  developed  and  applied  to  study  inflammation 
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(23-33),  their  limited  scope  and  focus  predominantly  on  qualita¬ 
tive  representations  of  inflammation  dynamics  generally  restrict 
their  ability  to  provide  accurate  interpretations  of  existing  data 
sets  and  generate  (semi)quantitative  hypotheses.  In  this  article,  we 
introduce  a  kinetic,  inherently  quantitative  computational  model 
of  inflammation  whose  parameters  were  derived  directly  from 
in  vitro  data.  We  validated  the  model  using  experimental  data  on 
acute  inflammation  and  then  extended  it  to  investigate  chronic  in¬ 
flammatory  scenarios. 

Our  modeling  results  suggest  that  the  modulation  of  two 
mechanistic  parameters,  macrophage  influx  and  efflux  rates,  may 
be  the  strongest  general  trigger  of  chronic  inflammation.  Moreover, 
we  predicted  that  local  concentration  changes  of  three  molecular 
mediators — IL-6,  TGF-(3,  and  PDGF — may  be  sensitive  and  ro¬ 
bust  indicators  of  the  ongoing  or  imminent  chronic  inflammatory 
response  under  various  scenarios.  Our  results  are  in  agreement 
with  existing  experimental  and  clinical  data  and  suggest  that  key 
dynamic  features  of  the  inflammatory  response  can  be  explained 
and  predicted  using  a  parsimonious  computational  modeling  ap¬ 
proach. 


P  or  IL-lp  activation  and  use  the  simplifying  assumption  that,  because 
inflammation  is  present,  any  latent  TGF-p  or  IL-ip  that  is  released  will  be 
in  its  activated  form. 

Estimation  of  production  rates 

Among  its  69  main  parameters,  our  computational  model  contains  26 
production  rates.  Using  published  data  from  in  vitro  experiments,  the 
production  rates  were  estimated  as  follows: 


where  k  is  the  production  rate  for  a  cytokine  or  growth  factor,  C  is  the 
concentration  of  the  cytokine  (growth  factor)  measured  in  the  supernatant 
from  the  cell  culture  after  incubation,  M  is  the  initial  concentration  of 
macrophage  suspension  in  the  cell  culture  experiment,  and  t  is  the  time  of 
incubation  of  macrophages  with  a  stimulant  (e.g.,  bacterial  load  or  endo¬ 
toxin).  In  cases  where  measurements  were  available  for  multiple  time 
points,  linear  regression  was  used  to  calculate  slopes,  which  were  subse¬ 
quently  divided  by  the  initial  macrophage  concentration  to  give  the  final 
production  rate  values.  As  an  illustration,  Fig.  2 A  and  2B  show  the  IL-lp, 
TNF-a,  and  IL-10  production  rate  estimation  for  pro-  and  anti-inflammatoiy 
macrophages.  To  approximate  the  chemotactic  migration  data  for  neutrophils 
and  monocytes  (38,  39),  we  used  a  combination  of  quadratic  and  linear 
functions  with  at  least  three  parameters  (Fig.  2C,  2D). 


Materials  and  Methods 

Computational  model  and  simulations 

Our  mathematical  model  represents  local  inflammation  in  a  wound.  The 
model  consists  of  15  ordinary  differential  equations  and  1  delay  differential 
equation  (DDE;  Table  I),  and  it  has  69  main  parameters  representing 
different  molecular  and  cellular  processes,  such  as  neutrophil  and  mac¬ 
rophage  phenotype  conversion,  phagocytosis  of  apoptotic  neutrophils,  and 
the  production  and  degradation  of  cytokines  and  growth  factors  (Supplemental 
Table  I;  note  that  the  19  chemotaxis  function  parameters  in  the  model  are 
regarded  as  a  separate  group  and  are  not  included  among  the  main 
parameters).  In  the  model,  we  included  the  cell  types,  cytokines,  and  growth 
factors  that  are  commonly  regarded  as  key  components  of  the  inflammatory 
response  (Fig.  1).  The  DDE  in  the  model  is  used  to  describe  pro-inflam¬ 
matory  macrophage  dynamics.  In  the  bloodstream,  macrophages  exist  in  the 
form  of  their  precursors  (i.e.,  monocytes),  which  migrate  at  rates  similar  to 
those  for  neutrophil  migration  (34)  and  then  differentiate  into  active  mac¬ 
rophages  at  the  site  of  inflammation  within  — ^12 — 48  h  (16).  Because  we  did 
not  explicitly  define  monocyte  kinetics  in  our  model,  we  represented 
monocyte  differentiation  by  incorporating  a  12-h  delay  in  the  argument  of 
the  chemotaxis  function  that  drives  pro-inflammatory  macrophage  migra¬ 
tion  in  the  model  (33)  (Table  I). 

The  time  courses  for  the  concentrations  of  these  cell  types  and  molecular 
species  constitute  the  output  of  our  model.  Its  main  input  is  the  initial 
platelet  concentration  (default  value:  2  X  108  platelets/ml),  which  repre¬ 
sents  an  injury  initiating  local  inflammation  (all  other  model  variables  are 
initially  set  to  0).  The  model  equations  were  solved  numerically  using  the 
MATLAB  solver  DDE23  with  default  tolerance  values.  Our  model  is  in¬ 
trinsically  quantitative,  because  all  of  its  variables  and  main  parameters 
(except  the  dimensionless  feedback  functions  discussed  below)  are  expressed 
in  absolute  units,  and  the  majority  of  the  parameter  values  were  derived  di¬ 
rectly  from  experimental  data.  We  performed  all  computations  in  MATLAB 
2012a  (MathWorks,  Natick,  MA).  Each  simulation  covered  a  20-d  time  in¬ 
terval  after  inflammation  initiation. 

In  the  model,  we  assume  that  the  inflammatory  response  is  initiated  only 
by  TGF-p  released  by  platelets  postinjury.  Although  the  inflammatory 
response  also  can  be  stimulated  by  other  locally  secreted  chemoattractants 
[such  as  PDGF,  TNF-a,  MIP-la,  and  CXCL8  (35-37)],  TGF-(3  has  been 
characterized  as  the  strongest  inflammatory  cell  chemoattractant  in  vitro 
(38,  39)  and  is  used  in  our  model  as  a  proxy  for  all  initiating  chemotactic 
signals.  In  the  model,  we  assume  that  anti-inflammatory  macrophages  are 
generated  by  phenotype  conversion  of  pro-inflammatory  macrophages 
triggered  by  phagocytosis  of  apoptotic  neutrophils,  which  is  suggested  by 
experimental  data  (40-42).  Yet  another  modeling  assumption  pertains  to 
the  activity  of  TGF-(3  and  IL-lfT  When  released  from  various  cellular 
sources,  TGF-(3  exists  in  a  latent  state  (43)  and  is  subsequently  activated  by 
the  presence  of  injury-  and  infection-related  enzymes,  as  well  as  the  acidic 
and  elevated-temperature  environment  of  the  inflamed  site  (44).  Analo¬ 
gously,  the  biological  activity  of  IL-lp  requires  its  maturation  into  an 
activated  form,  which  occurs  during  or  after  secretion  through  the  cell 
membrane  and  is  preceded  by  inflammasome-mediated  cleavage  of  the 
translated  IL-ip  (45,  46).  We  do  not  explicitly  define  the  kinetics  of  TGF- 


Estimation  of  degradation  rates 

The  model  contains  12  degradation/removal  rates.  Half-life  estimates  for  the 
catabolic  breakdown  of  cytokines  and  growth  factors  were  obtained  from 
the  literature.  Assuming  a  first-order  decay  of  the  molecular  species,  the 
degradation  rates  were  calculated  as  follows: 


k  = 


0.693 

tl/2 


(2) 


where  tm  is  the  half-life  of  a  cytokine  or  growth  factor,  and  k  is  the 
corresponding  degradation  rate. 


Estimation  of  the  parameters  describing  phagocytosis  of 
apoptotic  neutrophils  by  pro-inflammatory  macrophages 

Pro-inflammatory  macrophages  undergo  a  change  in  phenotype  upon 
phagocytizing  apoptotic  neutrophils  and  cellular  debris  (40-42).  In  our 
model,  the  phagocytosis  rate  of  apoptotic  neutrophils  and  phenotype 
conversion  rate  for  macrophages  were  determined  by  approximating  ex¬ 
perimental  phagocytosis  data  using  hyperbolic  functions.  For  example,  the 
data  by  Newman  et  al.  (47)  characterize  the  number  of  apoptotic  neu¬ 
trophils  ingested  after  4  h  of  incubation  with  a  fixed  initial  concentration  of 
macrophages  for  different  initial  concentrations  of  apoptotic  neutrophils. 
These  data  were  converted  into  rates  of  apoptotic  neutrophil  ingestion  per 
pro-inflammatory  macrophage  per  hour  for  different  initial  apoptotic 
neutrophil  concentrations  (Supplemental  Fig.  1A)  and  were  then  fitted 
with  the  following  equation: 


9  ( Nap0p )  — 


^1  ingest  Napop 
klingest  T  N apop 


(3) 


where  Napop  is  the  concentration  of  apoptotic  neutrophils,  and  k\ingest  and 
kiingest  are  the  model  parameters  that  were  estimated  as  a  result  of  the 
fitting  (Supplemental  Table  I).  The  fitting  was  performed  using  MAT- 
LAB’ s  curve-fitting  toolbox. 

An  expression  similar  to  Eq.  3  was  used  to  describe  the  rate  of  phenotype 
conversion  of  pro-inflammatory  into  anti-inflammatory  macrophages, 
which  occurs  as  a  result  of  apoptotic  neutrophil  phagocytosis  by  pro- 
inflammatory  macrophages.  However,  phagocytosis  is  not  a  “one-on-one 
process”  (i.e.,  a  single  macrophage  can  phagocytize  up  to  three  apoptotic 
neutrophils)  (41,  48,  49).  Using  the  data  reported  by  Newman  et  al.  (47), 
we  calculated  the  fraction  of  macrophages  ingesting  apoptotic  neutrophils 
per  hour  for  different  initial  concentrations  of  apoptotic  neutrophils 
(Supplemental  Fig.  IB).  These  data  were  then  fitted  with  the  following 
equation: 


9  (Napop) 


klingest  Nap0p 

klingest  T  Napop 


(4) 


where  kungest  and  k2ingest  are  the  model  parameters  that  were  estimated  as 
a  result  of  the  fitting  (Supplemental  Table  I). 
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CHRONIC  INFLAMMATION  MECHANISMS  AND  MOLECULAR  INDICATORS 


Estimation  of  cytokine  feedback  functions 

Pro-  and  anti-inflammatory  macrophages,  neutrophils,  and  platelets  secrete 
cytokines  and  growth  factors  that  provide  regulatory  feedback  for  the  in¬ 
flammatory  response  by  upregulating  (positive  feedback)  or  downregulating 
(negative  feedback)  the  production  of  other  cytokines.  (We  generally  regard 
these  interactions  as  feedbacks  because  they  represent  the  self-modulating 
effects  of  extracellular  mediators  on  the  intracellular  machinery  responsible 
for  the  production  of  such  mediators.)  To  reflect  these  effects  in  the  model, 
we  introduced  10  dimensionless  feedback  functions  denoted  /],  f2, 
(Supplemental  Table  I)  that  represent  fractional  increases  or  decreases 
(induced  by  a  particular  cytokine  or  growth  factor)  in  the  production  rates 
of  other  cytokines  and  growth  factors  for  pro-inflammatory  macrophages. 
The  parameters  of  these  functions  were  estimated  by  fitting  the  functions 
to  experimental  data  using  the  curve-fitting  toolbox  in  MATLAB 
(Supplemental  Fig.  1C). 

Sensitivity  analysis 

We  calculated  logarithmic  (i.e.,  relative)  local  sensitivities,  s^t),  at  time 
moments  t,  according  to  the  standard  definition  (see,  e.g.,  Ref.  50): 

Sij(t)  =  dlogXi^/dlogpj  =  (dXi/Xi)/(d pj/pj)  (5) 

where  Xft)  is  the  model’s  ith  variable  and  pj  is  the  model’s  j'th  parameter 
(of  the  model’s  69  main  parameters).  By  definition,  each  of  these  sensi¬ 
tivities  reflects  the  magnitude  of  the  relative  change  in  a  model’s  output 
variable  induced  by  a  local  (i.e.,  small)  relative  change  in  a  model’s  pa¬ 
rameter.  To  obtain  numeric  approximations  of  the  derivatives  in  Eq.  5,  each 
parameter  was  individually  perturbed  by  ±1%  of  its  value,  and  the  de¬ 
rivative  was  approximated  using  the  second-order  central  finite  difference 
formula.  We  performed  local  sensitivity  analysis  in  the  vicinity  of  the 
default  parameter  set,  as  well  as  for  10,000  random  parameter  sets  in 
which  individual  parameters  were  sampled  independently  from  intervals, 
permitting  up  to  2-fold  deviations  (up  or  down)  from  the  corresponding 
default  values.  This  random  sampling  was  intended  to  represent  possible 
natural  variations  in  the  molecular  environment  of  the  inflamed  site  for 
different  subjects  under  different  inflammatory  scenarios.  To  generate  the 
random  parameter  sets,  we  used  Latin  hypercube  sampling,  as  previously 
described  (51)  (we  used  the  MATLAB  function  LHSDESIGN).  In  all 
analyses,  we  calculated  local  sensitivities  for  each  of  the  21  evenly  spaced 
time  points  that  discretize  the  total  20-d  simulation  interval  into  1-d  sub¬ 
intervals  (i.e.,  day  0,  day  1,  and  so  forth).  To  compare  and  rank  sensitiv¬ 
ities,  we  used  their  absolute  values. 

Results 

Model  captures  essential  kinetic  features  of  acute  inflammation 

We  tested  our  model’s  ability  to  capture  typical  features  of  the 
time  course  of  acute  inflammation  by  comparing  modeling  pre¬ 
dictions  with  published  experimental  data.  The  model  includes 
descriptions  for  key  mechanisms  involved  in  the  inflammatory 
response  (Fig.  1),  with  the  parameters  quantifying  these  mecha¬ 
nisms  derived  from  acute  inflammatory  response  data  (Fig.  2).  The 
model  predicted  the  occurrence  of  a  neutrophil  peak  at  ~1  d  and 
a  macrophage  peak  at  ~2  d  after  inflammation  initiation  (Fig.  3 A, 
3B).  These  predictions  are  in  accord  with  experimental  measure¬ 
ments  showing  that  a  single,  well-pronounced  peak  typically 
occurs  at  1  —  3  d  after  inflammation  initiation  for  neutrophils  and 
at  2—4  d  for  macrophages  (3,  52).  Furthermore,  it  is  expected  that 
the  concentrations  of  the  inflammatory  cells  and  cytokines  return 
to  their  basal  level  within  2—3  wk  after  the  onset  of  inflammation 
(3,  4,  52).  Our  model  predicted  the  return  of  all  variables  in  the 
model  to  their  default  values  within  20  d  after  inflammation  ini¬ 
tiation  (Fig.  3A-3C).  Moreover,  our  model  predicted  the  kinetics 
of  individual  neutrophil  and  macrophage  phenotypes  (Supplemental 
Fig.  2). 

To  further  validate  our  computational  model,  we  compared  its 
predictions  with  the  dynamics  of  acute  inflammatory  response  in 
a  number  of  animal  and  human  in  vivo  models  (those  datasets  were 
not  used  to  estimate  the  model  parameters).  The  computational 
model  successfully  predicted  key  features  (such  as  the  overall  curve 
shape,  peak  time,  and  resolution  time)  for  neutrophil  and  macro¬ 


phage  time  courses,  as  observed  in  human  wounds  (52),  and  for 
peritoneal  infection  in  mice  (16)  (Fig.  3A,  3B).  We  validated  TGF- 
(3  concentration  predictions  using  experimental  data  from  rats  with 
excisional  wound  injury  (53)  (Fig.  3C).  Most  of  the  cytokines  and 
chemokines  in  our  model  showed  good  qualitative  or  semiquan- 
titative  agreement  with  available  inflammation  data.  As  an  illus¬ 
tration,  we  show  modeling  predictions  and  their  experimental 
validation  (with  experimental  data  taken  from  Ref.  54)  for  two 
representative  cytokines  with  pro-inflammatory  (IL-1J3)  and  anti¬ 
inflammatory  (IL-10)  properties  in  Fig.  3D  and  3E.  Interestingly, 
for  two  of  the  cytokines  [IL-12  (Fig.  3F)  and  IL-6  (data  not 
shown)],  our  model  gave  quantitatively  accurate  kinetic  predictions. 
Opportunities  for  a  more  extensive  quantitative  validation  were 
limited  as  the  result  of  interassay  and  interlaboratory  variations  in 
the  measured  concentrations  of  inflammatory  cells  and  molecular 
mediators. 

Modeling  predicts  macrophage  influx  and  efflux  modulation  to 
be  the  strongest  mechanistic  trigger  of  chronic  inflammation 

Our  analysis  of  chronic  inflammation  was  based  on  the  general 
assumption  that  delayed  resolution  of  inflammation  results  from 
dysregulation  of  certain  key  mechanisms  that  are  represented  by  the 
main  parameters  in  our  model.  We  used  this  assumption  to  extend 
our  acute  inflammation  model  to  represent  chronic  inflammatory 
scenarios.  This  was  achieved  by  performing  an  extensive  sensitivity 
analysis  to  identify  the  parameters  that  may  display  a  strong 
functional  association  with  chronic  inflammation. 

Experimental  data  suggest  that  the  pro-inflammatory  cytokines 
TNF-a,  IL-1(3,  and  IL-6  are  consistently  and  considerably  in¬ 
creased  in  chronic  inflammation  induced  by  infection  or  disease, 
such  as  diabetes  (12,  14-17).  We  initially  used  local  sensitivity 
analysis  in  the  vicinity  of  the  default  parameter  set  (see  Materials 
and  Methods)  to  identify  the  parameters  whose  change  induced 
the  largest  variation  in  the  concentrations  of  TNF-a,  IL-1(3,  and 
IL-6  (analyzed  independently).  We  found  that,  for  each  of  the 
three  cytokines,  the  largest  relative  changes  (for  the  majority  of 
the  time  points  considered  in  the  sensitivity  analysis)  were  in¬ 
duced  by  modulating  the  same  two  parameters  of  the  total  69 
main  parameters.  These  two  parameters  represented  the  pro- 
inflammatory  macrophage  influx  rate  and  the  rate  of  efflux  of  both 
pro-  and  anti-inflammatory  macrophages  (designated  as  kM  in  and 
kd  m  in  the  model,  respectively;  Table  I,  Supplemental  Table  I). 
For  earlier  time  points  in  the  simulation  (~2-10  d  after  inflam¬ 
mation  initiation),  macrophage  influx  rate  had  the  highest  impact 
on  the  cytokine  concentrations  (Fig.  4A),  whereas  macrophage 
efflux  rate  had  the  second  highest  influence  on  those  outputs  (data 
not  shown).  At  later  time  points  (~7-20  d  after  inflammation 
initiation),  the  macrophage  efflux  rate  emerged  as  the  strongest 
modulator  of  the  levels  of  all  three  cytokines  (Fig.  4A),  and  the 
macrophage  influx  rate  was  the  second  strongest  (data  not  shown). 
Overall,  these  two  parameters  represented  the  dominant  factors 
leading  to  the  overproduction  of  TNF-a,  IL-1(3,  and  IL-6  in  our 
model. 

We  tested  the  robustness  of  these  results  by  performing  local 
sensitivity  analysis  in  the  vicinities  of  10,000  randomly  selected 
parameter  sets.  Specifically,  we  identified  the  parameters  that  in¬ 
duced  the  largest  variations  in  the  concentrations  of  the  three 
cytokines  (i.e.,  TNF-a,  IL-1(3,  and  IL-6)  across  10,000  simulations 
with  random  parameter  sets  (see  Materials  and  Methods).  We 
verified  that,  for  >70%  of  the  simulations,  the  sensitivities  for  the 
two  parameters  (i.e.,  macrophage  influx  and  efflux  rates)  were  the 
highest  for  the  majority  of  the  considered  simulation  time  points. 
(Specifically,  the  macrophage  influx  rate  was  the  most  sensitive 
parameter  at  days  4—9  after  inflammation  initiation  for  TNF-a;  at 
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FIGURE  1.  Inflammation  model  components  and 
interactions.  Platelets  initiate  the  wound-healing  pro¬ 
cess  by  supporting  blood  clotting  in  the  wound.  They 
also  secrete  growth  factors,  such  as  TGF-(3  and  PDGF, 
which  promote  the  migration  of  neutrophils  and  mono¬ 
cytes  from  surrounding  blood  vessels  (43,  86).  Monocytes 
mature  into  pro-inflammatoiy  macrophages  at  the  wound 
site.  Neutrophils  ingest  cellular  debris  and  the  invading 
bacteria  and,  as  a  result,  become  apoptotic  (1).  Pro-in¬ 
flammatory  macrophages  then  phagocytize  the  apoptotic 
neutrophils  (and  cellular  debris)  and  change  their  pheno¬ 
type  to  anti-inflammatory  (41).  Both  pro-  and  anti¬ 
inflammatory  macrophages  secrete  cytokines  and  growth 
factors  at  different  rates,  and  these  molecular  mediators 
provide  positive  and/or  negative  feedback  to  the  ongoing 
inflammation  (6,  87).  Anti-inflammatoiy  macrophages  are 
cleared  by  migration  to  lymphatic  vessels  (88). 


Cells 

tissue  0  Platelet 

Active  neutrophil 
Apoptotic  neutrophil 
Pro-inflammatory  macrophage 
Anti-inflammatory  macrophage 


Cytokines  and  growth  factors 
TNF-a,  IL-1p,  IL-6,  IL-10,  IL-12, 
CXCL8,  MIP-la,  MIP-2,  IP-10, 
TGF-P,  and  PDGF 


. *>  Chemotaxis/migration 

- ►  Secretion/activation 

- ►  Differentiation 
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wound  site 


days  2,  6,  and  7  for  IL- 1  (3 ;  and  at  days  2-10  for  IL-6.  The 
macrophage  efflux  rate  was  the  most  sensitive  parameter  at  days 
10—20  after  inflammation  initiation  for  TNF-a;  at  days  8—20  for 
IL-1(3;  and  at  days  11—20  for  IL-6.)  Representative  simulation 
results  for  3  d  (days  4,  10,  and  16)  are  shown  in  Fig.  4B.  For  each 
of  the  10,000  randomly  selected  parameter  sets,  just  like  for  the 
default  parameter  set,  macrophage  influx  rate  had  the  largest  in¬ 
fluence  on  all  three  cytokines  at  earlier  simulation  time  points, 
whereas  macrophage  efflux  rate  dominated  at  later  times.  This 
result  is  consistent  with  the  dominant  role  of  macrophage  influx 
during  inflammation  initiation  and  the  prevalent  role  of  macro¬ 
phage  efflux  from  the  inflamed  site  during  inflammation  resolu¬ 
tion.  Notably,  the  two  parameters  induced  changes  in  opposite 
directions:  an  increase  in  the  concentrations  of  the  three  cytokines 


[TGF-P]  (pg/ml)  [TGF-p]  (pg/ml) 


FIGURE  2.  Estimation  of  cytokine  production  rates  and  chemotaxis 
functions  from  experimental  data.  Experimental  data  on  the  production  of 
IL-ip  (♦),  TNF-a  (■),  and  IL-10  (a)  were  obtained  for  cultured  human 
pro-inflammatory  (A)  and  anti-inflammatory  (B)  monocytes  (40)  and  ap¬ 
proximated  using  linear  regression  analysis  (solid  lines).  The  slopes  of  the 
regression  lines  represent  estimated  cytokine  production  rates.  Experi¬ 
mental  data  on  neutrophil  (•;  C)  and  monocyte  (O;  D)  chemotaxis  in¬ 
duced  by  TGF-p  were  obtained  from  Brandes  et  al.  (39)  and  Wahl  et  al. 
(38).  The  migration  rate  functions  for  both  neutrophils  and  monocytes 
were  obtained  by  fitting  the  experimental  data  using  a  combination  of 
quadratic  and  linear  functions  (solid  lines). 


could  be  caused  by  increasing  macrophage  influx  rate  or  by  de¬ 
creasing  macrophage  efflux  rate.  Taken  together,  our  results  sug¬ 
gest  that  macrophage  flux  modulation  may  be  a  principal  cause  of 
a  delay  in  inflammation  resolution.  This  general  prediction  is 
consistent  with  studies  showing  that  a  decrease  in  macrophage 
efflux  from  the  wound,  resulting  from  disrupted  synthesis  of  the 
integrin  Mac-1,  delayed  wound  healing  in  mice  (55,  56). 

Our  local  sensitivity  analysis  results  suggested  that  compara¬ 
tively  large  changes  in  macrophage  influx  and  efflux  rates  can 
induce  a  considerable  and  possibly  sustained  increase  in  the  levels 
of  inflammatory  cells,  which  may  be  indicative  of  chronic  in¬ 
flammation.  We  tested  this  prediction  by  increasing  the  default 
value  of  macrophage  influx  rate  by  5-fold,  which  led  to  a  noticeable 
increase  in  the  total  neutrophil  concentration  and  a  considerable 
increase  in  the  total  macrophage  concentration  (Fig.  5).  In  com¬ 
parison  with  an  acute  inflammatory  scenario,  macrophage  reso¬ 
lution  was  delayed  significantly,  as  is  expected  for  chronic 
inflammation.  We  validated  these  computational  predictions  by 
comparisons  with  experimental  data  from  an  in  vivo  mouse  model 
for  both  acute  and  delayed  inflammatory  response  induced  by  the 
endotoxin  zymosan  given  i.p.  (16)  (Fig.  5).  [Such  a  comparison 
was  justified  by  the  assumption  that  a  higher  dose  of  endotoxin 
administration  leads  to  an  increased  macrophage  influx  into  the 
inflamed  area  (14,  57)].  In  summary,  these  results  demonstrate  that 
the  processes  of  macrophage  influx  and  efflux  may  be  critical  to 
a  well-balanced  inflammatory  response  and  that  their  dysregulat- 
ion  is  a  strong  trigger  of  chronic  inflammation. 

Modeling  identifies  IL-6,  TGF-/3,  and  PDGF  as  sensitive  local 
indicators  of  chronic  inflammation 

Having  established  that  macrophage  influx  and  efflux  rates  may 
play  a  defining  role  in  the  genesis  of  chronic  inflammation,  we 
wanted  to  determine  which  of  the  1 1  model  variables  representing 
inflammatory  mediators  would  be  most  sensitive  to  the  modulation 
of  macrophage  flux  rates.  We  first  analyzed  the  logarithmic  sen¬ 
sitivities  for  these  1 1  variables  (one  at  a  time)  with  respect  to  the 
two  parameters  (i.e.,  macrophage  influx  and  efflux  rates)  for  the 
default  parameter  set.  This  analysis  showed  that,  at  each  of  the 
considered  simulation  time  points  with  the  exception  of  day  1  after 
inflammation  initiation,  the  concentrations  of  IL-6,  TGF-(3,  and 
PDGF  were  characterized  by  the  highest,  second-highest,  and 
third-highest  sensitivities,  respectively,  for  both  of  the  perturbed 
macrophage  flux  parameters.  This  finding  suggested  that  these  three 
variables  could  be  sensitive  indicators  of  chronic  inflammatory 
conditions  induced  by  variations  in  the  macrophage  flux  rates. 
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CHRONIC  INFLAMMATION  MECHANISMS  AND  MOLECULAR  INDICATORS 


FIGURE  3.  Modeling  predictions  capture  experi¬ 
mentally  detected  kinetics  of  acute  inflammation.  Solid 
lines  show  modeling  predictions;  symbols  show  ex¬ 
perimental  data.  Brackets  designate  concentration.  (A- 
C)  Experimental  data  for  mouse  infection  (■)  (16),  rat 
wounds  (♦)  (53),  and  human  wounds  (•)  (52).  (D  and 
E)  Normalized  values  of  IL- 1  (3  and  IL-10  concen¬ 
trations  in  vitro.  Experimental  data  from  cultured  hu¬ 
man  macrophages  are  shown  (a)  (54).  For  proper 
comparisons  between  model  predictions  and  experi¬ 
mental  data,  normalization  was  necessary  because  of 
differences  in  reporting  units  between  experimental 
data  and  model  outputs  and  was  achieved  by  division  of 
all  of  the  predicted  or  experimentally  obtained  values  of 
a  particular  variable  by  the  maximum  value  of  that 
variable.  (F)  IL-12  concentration  dynamics  in  vitro. 
Experimental  data  were  taken  from  Kum  et  al.  (54). 


To  test  the  robustness  of  these  indicators,  we  analyzed  the 
logarithmic  sensitivities  for  the  10,000  randomly  selected  pa¬ 
rameter  sets.  For  each  of  the  considered  21  simulation  time  points 
and  each  of  the  11  variables,  we  calculated  the  fractions  of  the 
random  parameter  sets  for  which  a  given  variable  demonstrated  the 
highest,  second-highest,  or  third-highest  sensitivity  with  respect  to 
macrophage  influx  or  efflux  rate  perturbations.  Figs.  6  and  7  show 
the  results  for  4  d  in  the  simulation  (days  2,  4,  10,  and  16  after 
inflammation  initiation).  These  fractions  were  then  averaged  over 
the  first  10  d  after  inflammation  initiation,  because  we  focused  on 
the  clinically  important,  early  phases  of  inflammation  develop¬ 
ment.  We  found  that  IL-6  concentration  ranked  as  the  most  sen¬ 
sitive  variable  for  an  average  of  84%  of  the  randomly  selected 
parameter  sets  when  the  perturbed  parameter  was  macrophage 
influx  rate  and  for  an  average  of  67%  of  the  randomly  selected 
parameter  sets  when  the  perturbed  parameter  was  macrophage 
efflux  rate.  TGF-(3  concentration  ranked  as  the  second  most  sen¬ 
sitive  variable  for  66  and  61%  (on  the  average)  of  the  parameter 
sets  when  the  perturbed  parameters  were  macrophage  influx  rate 
and  macrophage  efflux  rate,  respectively.  Likewise,  for  these  two 
flux  parameters,  PDGF  concentration  ranked  as  the  third  most 
sensitive  variable  for  68  and  50%  of  the  parameter  sets,  respec¬ 
tively. 

These  results  demonstrate  that,  among  all  of  the  inflammatory 
mediators  represented  in  the  model,  IL-6,  TGF-(3,  and  PDGF  had 
the  largest  relative  changes  induced  by  chronic  inflammation 
triggers.  Although  the  levels  of  TNF-a  and  IL-1J3  [often  reported 
in  the  literature  as  chronic  inflammation  indicators  (12,  14-17)] 
also  can  be  increased  by  macrophage  flux  modulation,  the  relative 
magnitude  of  those  increases  in  our  modeling  predictions 
appeared  to  be  generally  smaller  than  those  for  IL-6,  TGF-J3,  and 
PDGF.  Our  findings  suggest  that  these  three  molecular  mediators 
(particularly  IL-6)  can  be  used  as  reliable  early  indicators  of 
chronic  inflammation  in  clinical  settings.  Interestingly,  this  con¬ 
clusion  is  consistent  with  a  recent  clinical  study  of  combat  wounds 
by  Hawksworth  et  al.  (19),  which  identified  local  IL-6  level  as  the 
best  predictor  of  dehiscence  in  traumatic  wounds. 

Model  robustness  analysis 

We  verified  the  robustness  of  the  constructed  computational  model 
in  the  vicinity  of  the  default  parameter  set  using  local  sensitivity 
analysis,  as  described  in  Materials  and  Methods.  Specifically,  we 
confirmed  that  none  of  the  model’s  16  output  variables  was  overly 


sensitive  [i.e.,  Syit)  >  10]  to  any  of  the  parameter  perturbations. 
Because  a  major  source  of  complexity  in  the  inflammatory  system 
is  the  functional  (feedback)  interactions  between  different  in¬ 
flammatory  mediators,  we  chose  to  limit  such  feedbacks  repre¬ 
sented  in  the  model  (Supplemental  Table  I)  to  only  those  whose 
elimination  (modeled  by  setting  the  corresponding  feedback 
functions  to  1)  significantly  (i.e.,  by  >5%)  affected  at  least  one  of 
the  model  outputs.  One  such  feedback,  present  in  an  earlier  ver¬ 
sion  of  the  model,  represented  the  downregulation  of  the  pro¬ 
duction  of  the  pro-inflammatory  chemokine  CXCL8  by  the  anti¬ 
inflammatory  cytokine  IL-10.  This  functional  connection  did  not 
meet  the  above  selection  criterion  and  was  excluded  from  the 
current  version  of  the  model.  Interestingly,  all  of  the  major  results 
reported  in  the  preceding  subsections  for  the  current  version  of  the 
model  also  held  for  the  version  in  which  that  feedback  interaction 
was  present. 

We  found  that  elimination  of  a  single  feedback  induced  a  >5% 
change  only  in  the  variable  that  the  feedback  directly  modulated 
(i.e.,  its  “affected”  variable;  see  Supplemental  Table  I).  However, 
it  is  conceivable  that  eliminating  all  feedbacks  associated  with  one 
or  more  “effector”  variables  (i.e.,  the  independent  variable  of  the 
corresponding  feedback  functions;  see  Supplemental  Table  I) 
could  considerably  influence  other  variables  in  addition  to  the 
feedbacks’  “affected”  variables.  To  investigate  this  possibility,  we 
performed  simulations  in  which  we  simultaneously  eliminated  all 
feedbacks  driven  by  the  anti-inflammatory  mediators  TGF-(3  and 
IL-10  (feedback  functions  /4,  /5,  /9,  /10  and  /j,  /2,  /3,  respectively; 
see  Supplemental  Table  I)  for  the  acute  and  chronic  inflammatory 
scenarios.  Indeed,  elimination  of  these  two  groups  of  feedbacks 
induced  a  >5%  change  in  the  neutrophil  and  macrophage  con¬ 
centrations  for  both  acute  and  chronic  inflammation  (Supplemental 
Fig.  3C,  3D),  demonstrating  the  importance  of  interactions 
among  feedbacks.  The  resulting  changes  for  chronic  inflamma¬ 
tion  were  comparatively  large  (up  to  ~50  and  ~100%  for  neu¬ 
trophils  and  macrophages,  respectively;  compare  dashed  and 
dash-dotted  lines  in  Supplemental  Fig.  3C,  3D)  and  exceeded 
those  for  acute  inflammation.  This  indicates  that  functional 
contributions  of  feedbacks  can  have  increased  significance  when 
the  inflammatory  process  is  pathologically  perturbed. 

To  investigate  the  effects  of  feedback  control  exerted  by 
pro-inflammatory  mediators,  we  performed  a  simulation  in  which 
TNF-a  (and  therefore  all  of  the  associated  feedbacks)  was 
“knocked  out”  (which  was  modeled  by  setting  all  of  the  TNF-a 


Downloaded  from  http://www.jimmunol.org/  at  UNIFORMED  SERVICES  UNIVERSITY  OF  THE  HEALTH  SCIENC  on  February  11,  2014 


The  Journal  of  Immunology 


1829 


Table  I.  Model  variables  and  equations 


Model  variables:  volume  concentrations  for  different  cell  types  and  molecular  species 


X act 

Active  neutrophils 

TGFp 

Transforming  growth  factor- P 

N 

1  v  apop 

Apoptotic  neutrophils 

PDGF 

Platelet-derived  growth  factor 

FI pro 

Pro-inflammatory  macrophages 

IL\p 

Interleukin- 1  p 

FI  anti 

Anti-inflammatory  macrophages 

IL6 

Interleukin- 6 

P 

Platelets 

FIIP  \  a 

Macrophage  inflammatory  protein- la 

CXCLS 

Chemokine  CYCL8 

FIIP2 

Macrophage  inflammatory  protein-2 

IL\2 

Interleukin- 12 

IP\0 

Interferon-y-induced  protein  10 

7Z10 

Interleukin- 10 

TNFa 

Tumor  necrosis  factor-a 

Model  equations 


dP  /  dt  =  -kd  pP; 

dNact  /  dt  =  kNJn  {fN  ( TGFJ3 )  +  fN  (PDGF)  +  f(CXCLK)}  -  kN_deactNact ; 


dN apopt  I  dt  —  kN  deactiy  act 


k  N 

JSJ  _  \ingest  apop 

-  k  +  N  pro  ’ 

Hugest  apop 


dM  /dt  =  kM  {fM (TGFfKt - 12))  +  fM (PDGF)  +  f(TNFa)  +  f(MIP\a)\ 


r  r  ,T  \ 

k  V 

Hugest  apop 
y  ^ Hugest  ^  apop  J 


M  pro  kd  M 


f  k  N  ^ 

Y  ^  Hugest  apop 

^ Hugest  ^  apop  J 


M 


pro  ’ 


dMauti  / dt  =  ((kUngestNapop )  /(k2ingest  +  Napop  ))Mpro  -  kd  MM anti ; 

dTNFa  /  dt  =  kTNFa actN act  +  kTNFa  profxfAf6M pro  +  kTNFa  antiManti  —  kd  mFaTNF(X ; 

dTGF/3  /  dt  =  kT  pP + kTGFp  proM pro  +  kTGFp  antiManti  —  kd  TGFpTGFfi ; 


dPDGF  /  dt  —  kPDGF  proM pro  +  kPDGF  antiM anti  —  kd  pdgfPDGF  ; 

dILX/3  /  dt  —  kILXp  actNact  +  ^ i /? _  pro /3 /5 /7 M pro  +  kILlp  antiManti  —  kd  ILlpIL\(3 ; 

dCXCLR  /  dt  =  kCXCL%  proM pro  +  kCXCLS  antiM anti  —  kd  CXCL 8  CXCL 8 ; 

dIL6 1 dt  =  kIL6  actNact  +  kIL6  prof2 (1  +  f9)Mpro  +  kIL6  antiManti  -  kd  IL6IL6 ; 

r//Zl  0/  dt  =  kIL10  pro  (1  +  /10  )M pro  +  kILl0  antiManti  —  kd  IL]()IL\  0 ; 

J/Ll  2  /  dt  —  kILl2  profsMpro  —  kd  iluIL\  2 ; 


dMIPla  I  dt  —  kMIPla  proM pro  +  kMIPXa  antiM anti  —  kd  MIPXaMIP\a ; 
dMIP2 1  dt  =  kM1P2  proMpro  +  kMIP2  antiM anti  —  kd  M1P2MIP2 ; 

J/Pl  0  /  dt  =  kIPXO  proMpro  +  kjPXO  antiM anti  —  kd  IPX0IPl  0 . 


All  cell  concentration  variables  in  our  model  are  expressed  in  units  of  cells/ml.  Variables  for  molecular  mediators  (cytokines, 
chemokines,  and  growth  factors)  are  expressed  in  units  of  ng/ml.  Time  is  expressed  in  hours.  Feedback  and  chemotaxis  functions 
for  different  cytokines  were  numbered  and  are  described  in  detail  in  Supplemental  Table  I. 


production  rates  to  0).  The  simulated  TNF-a  knockout  resul¬ 
ted  in  decreased  concentration  peaks  for  neutrophils  and  macro¬ 
phages,  which  is  consistent  with  the  experimentally  observed 
decrease  in  the  leukocyte  levels  in  the  wounds  of  the  TNF-a 
receptor  p55-knockout  mice  (58)  (Supplemental  Fig.  3A,  3B). 
Taken  together,  our  findings  support  the  notion  that  feedback 
regulation  is  important  for  fine-tuning  the  inflammatory  re¬ 
sponse. 

The  first  working  version  of  our  model  comprised  only  11 
variables  and  28  parameters  and  was  much  simpler  than  the  current 
version  (which  is  discussed  in  this  article).  It  did  not  contain  the 
positive-  and  negative-feedback  interactions  between  the  cytokines 
shown  in  Supplemental  Table  I.  Moreover,  it  described  the  che- 
motactic  effects  of  only  two  (TGF-J3  and  PDGF)  of  the  five  (TGF- 
P,  PDGF,  CXCL8,  TNF-a,  and  MIP-la)  molecular  mediators 
whose  chemotactic  effects  are  reflected  in  the  current  version  of 
the  model.  Yet,  perhaps  surprisingly,  even  that  simpler  model  gave 


kinetic  predictions  that  were  in  good  agreement  with  experimen¬ 
tally  detected  acute  inflammatory  behavior  (data  not  shown).  This 
finding  suggests  that,  although  the  current  version  of  the  model  is 
a  simplified  representation  of  the  complexity  of  the  inflammatory 
process  in  vivo,  our  omission  of  many  of  its  functional  details  did 
not  seem  to  have  precluded  the  generation  of  reasonable  pre¬ 
dictions.  Furthermore,  local  sensitivity  analysis  of  our  original 
model  predicted  the  same  two  mechanisms  crucial  for  triggering 
chronic  inflammation  as  the  more  comprehensive,  current  version 
of  the  model.  The  most  sensitive  indicator  of  chronic  inflamma¬ 
tion  predicted  by  the  original,  simpler  model  was  IL-6,  whereas 
TNF-a  and  TGF-(3  appeared  among  the  three  most  sensitive 
indicators.  These  results  suggest  that  key  properties  of  the  in¬ 
flammatory  process  might  be  largely  insensitive  to  the  details  of 
many  complex  interactions  occurring  in  the  system  and  can  be 
understood  by  considering  only  the  essential  molecular  and  cel¬ 
lular  mechanisms. 
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A 

■  Macrophage  efflux  rate  (PI 0) 

■  Other  parameters 


FIGURE  4.  Sensitivity  analysis  identifies  strong 
triggers  of  chronic  inflammation.  (A)  Parameters  in¬ 
ducing  the  largest  changes  in  TNF-a,  IL-ip,  and  IL-6 
concentrations  identified  by  local  sensitivity  analysis  in 
the  vicinity  of  the  default  parameter  set  (see  Materials 
and  Methods).  (B)  Summary  of  local  sensitivity  anal¬ 
ysis  results  for  10,000  randomly  selected  parameter  sets 
(see  Materials  and  Methods )  for  three  representative 
days  (4,  10,  and  16).  Shown  is  the  percentage  of  the 
10,000  simulations  for  which  the  macrophage  influx 
(blue)  and  efflux  (red)  parameters  [identified  in  the  anal¬ 
ysis  in  (A)]  induced  the  largest  changes  in  TNF-a, 
IL-10,  and  IL-6. 
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Discussion 

Chronic  inflammation  is  becoming  increasingly  recognized  as  a 
major  contributor  to  several  known  pathologies,  such  as  cancer, 
rheumatoid  arthritis,  and  diabetic  ulcers,  among  many  others  (10, 
11).  Yet,  mechanistic  causative  factors  of  chronic  inflammation,  as 
well  as  its  reliable  molecular  indicators,  have  not  been  sufficiently 
characterized.  In  this  study,  we  used  computational  modeling 
approaches  to  gain  mechanistic  insights  into  chronic  inflamma¬ 
tion.  Our  differential  equation-based  model  successfully  predicted 
key  kinetic  features  of  the  local  acute  and  chronic  inflammatory 
response,  as  was  verified  by  direct  comparisons  of  modeling 
predictions  with  experimental  data.  Using  sensitivity  analysis,  we 
predicted  that  the  levels  of  the  pro-inflammatory  cytokines  IL-6, 
IL-10,  and  TNF-a  are  affected  the  most  by  the  same  two  model 


FIGURE  5.  Modeling  predictions  capture  experimentally  observed 
chronic  inflammatory  kinetics.  Modeling  predictions  and  experimental 
validation  for  acute  inflammation,  taken  from  Fig.  3  (■),  are  shown  here 
for  comparison  purposes.  In  Fig.  3A-D,  the  computational  predictions  and 
experimental  data  were  normalized  (separately)  to  the  maximum  value. 
(Note  that  normalization  in  this  figure  is  different  from  that  for  Fig.  3 A  and 
3B,  because  the  maximum  values  in  the  data  sets  were  different.)  Solid 
(acute  inflammation)  and  dashed  (chronic  inflammation)  lines  show  model 
predictions.  ■  show  experimental  data  for  acute  inflammation;  #  show 
experimental  data  for  chronic  inflammation.  Chronic  inflammation  in 
the  computational  model  was  induced  by  increasing  macrophage  influx 
rate  (P7)  5-fold  from  its  default  value.  Chronic  inflammation  with  delayed 
healing  in  mice  was  obtained  by  a  100-fold  higher  endotoxin  dose 
given  i.p.  (16).  Shown  are  the  total  neutrophil  (A)  and  macrophage  (B) 
concentrations. 


parameters:  the  rates  of  macrophage  influx  into  and  efflux  out  of 
the  inflamed  site.  Because  the  levels  of  these  cytokines  are  typi¬ 
cally  elevated  in  chronic  inflammation,  our  results  suggest  that  the 
two  macrophage  flux  parameters  may  represent  the  mechanisms 
whose  malfunction  often  delays  normal  resolution  of  inflamma¬ 
tion.  Furthermore,  using  our  computational  model,  we  established 
that  different  model  variables  respond  differently  to  macrophage 
flux  modulation.  Of  the  11  model  variables  representing  inflam¬ 
matory  mediators,  the  concentrations  of  IL-6,  TGF-(3,  and  PDGF 
(in  order  of  decreasing  effect  magnitude)  were  most  affected  by 
macrophage  influx  and  efflux  rate  modulation.  This  finding  indi¬ 
cates  the  possibility  of  using  these  three  proteins,  which  can  be 
measured  in  the  local  environment  at  the  site  of  inflammation,  as 
reliable  predictive  and/or  diagnostic  indicators  of  the  development 
of  chronic  inflammation. 

Macrophages  are  regarded  as  the  most  multifunctional  immune 
cells,  influencing  nearly  all  phases  of  the  inflammatory  response. 
The  possibility  of  their  targeted  manipulation  has  naturally  been  the 
subject  of  extensive  investigations  (59-61).  Because  macrophages 
are  the  major  producers  of  cytokines  with  pro-  or  anti-inflammatory 
properties  (41,  59,  61),  any  significant  disruption  in  macrophage 
behavior  may  be  expected  to  alter  the  balance  of  cytokines  in  the 
inflammatory  milieu.  Elevated  levels  of  macrophage  infiltration  can 
lead  to  wound  fibrosis  (59).  At  the  same  time,  severe  early  inhibi¬ 
tion  of  macrophage  signaling  (as  well  as  macrophage  ablation 
during  later  stages  of  inflammation)  results  in  delayed  healing  in 
mouse  wound  models  (62,  63),  suggesting  that  initial  macrophage 
availability  (and,  therefore,  their  sufficient  influx)  is  necessary  for 
normal  inflammatory  response.  Thus,  our  modeling  prediction  re¬ 
garding  macrophage  flux  modulation  as  possibly  the  strongest 
overall  inducer  of  chronic  inflammation  can  perhaps  be  readily 
rationalized.  Nevertheless,  it  is  noteworthy  that,  in  our  model, 
macrophage  flux  rates  have  a  larger  impact  on  cytokine  production 
than  any  of  the  cytokines’  individual  production  and  degradation 
rates.  This  effect  might  be  (at  least  partially)  due  to  the  functional 
redundancy  exhibited  by  some  cytokines,  so  that  a  change  in  the 
production  or  degradation  rate  of  one  cytokine  can  be  compen¬ 
sated  for  by  the  action  of  the  others.  For  example,  TNF-a  and 
MIP-la  act  as  chemoattractants  for  macrophages  and  perform 
positive  regulation  of  the  production  of  other  pro-inflammatory 
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FIGURE  6.  Sensitivity  analysis  identifies  IL-6,  TGF-(3, 
and  PDGF  as  the  most  sensitive  indicators  of  chronic  in¬ 
flammation  induced  by  macrophage  influx  manipulation. 
Shown  are  the  percentages  of  the  total  10,000  random 
parameter  sets  for  which  different  molecular  species  in  the 
model  demonstrated  the  largest  (solid  bars),  second  largest 
(striped  bars),  and  third  largest  (stippled  bars)  relative 
change  induced  by  modulation  of  macrophage  influx. 
Model  predictions  are  shown  for  day  2  (A),  day  4  (B),  day 
10  (C),  and  day  16  (D)  after  inflammation  induction. 
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cytokines  (Supplemental  Table  I).  In  contrast,  because  of  a  lack 
of  compensatory  mechanisms,  a  disruption  in  the  migration  of 
pro-  and/or  anti-inflammatory  macrophages  has  a  direct  and 
global  impact  on  cytokine  production  and  downstream  pro¬ 
cesses. 

A  need  for  guided  therapy  has  been  well  documented  in  the 
treatment  of  chronic  wounds  resulting  from  military  trauma  (19, 
64).  A  treating  clinician’s  decision  to  surgically  close  a  wound  is 
commonly  based  on  subjective  and  somewhat  arbitrary  measures, 
such  as  the  patient’s  general  condition,  the  appearance  of  bones 
and  soft  tissues,  and  visual  signs  of  infection  or  necrosis  (65), 
which  are  poor  indicators  of  the  local  inflammatory  condition 
(66).  This  often  has  led  to  wound  dehiscence  and  nonhealing 
wounds.  Our  work  was  largely  motivated  by  the  need  to  find 
molecular  indicators  that  are  characteristic  of  the  local  wound 
environment  and  are  less  susceptible  to  systemic  influences 
and,  therefore,  better  reflect  ongoing  local  inflammation. 

Our  modeling  analysis  showed  that  three  molecular  mediators 
(IL-6,  TGF-P,  and  PDGF)  may  be  expected  to  demonstrate  con¬ 
sistently  large  relative  changes  in  their  local  concentrations  in¬ 


duced  by  independent  variations  in  macrophage  influx  or  efflux 
rates  (Figs.  6,  7).  In  our  analysis,  the  IL-6  level  ranked  as  the  most 
sensitive  indicator  of  local  changes  in  both  macrophage  flux 
parameters.  This  is  consistent  with  a  recent  clinical  study  by 
Hawksworth  et  al.  (19),  who  identified  IL-6  as  an  informative 
indicator  of  delayed  wound  healing.  However,  the  comparatively 
small  subject  group  and  limited  temporal  resolution  of  this  study 
indicate  the  need  for  additional  evidence.  Our  modeling  analysis, 
which  involved  an  extensive  parameter  randomization,  com¬ 
plements  these  clinical  results  and  suggests  that  IL-6  may  be 
considered  an  informative  molecular  indicator  of  chronic  in¬ 
flammation  for  a  wide  range  of  scenarios.  Furthermore,  the 
model-predicted  high  sensitivity  of  TGF-J3  to  macrophage  flux 
modulation  correlates  with  its  known  diverse  and  numerous 
roles  in  the  inflammation  process  (38,  39,  67-69).  Although  the 
clinical  study  by  Hawksworth  et  al.  (19)  did  not  report  results 
for  TGF-J3  and  PDGF,  our  computational  findings  may  provide 
motivation  for  experimental  and  clinical  testing  of  these  growth 
factors  as  additional  diagnostic  or  predictive  indicators  of  chronic 
inflammation. 
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FIGURE  7.  Sensitivity  analysis  identifies  IL-6,  TGF-(3, 
and  PDGF  as  the  most  sensitive  indicators  of  chronic  in¬ 
flammation  induced  by  macrophage  efflux  manipulation. 
Shown  are  the  percentages  of  the  total  10,000  randomized 
parameter  sets  for  which  different  molecular  species  in  the 
model  demonstrated  the  largest  (solid  bars),  second  largest 
(striped  bars),  and  third  largest  (stippled  bars)  relative 
change  induced  by  modulation  of  macrophage  efflux. 
Model  predictions  are  shown  for  day  2  (A),  day  4  (B),  day 
10  (C),  and  day  16  (D)  after  inflammation  induction. 
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The  importance  of  macrophage  flux  control  elucidated  in  this 
study  is  consistent  with  the  existing  view  of  macrophage  dynamics 
(i.e.,  monocyte  migration,  macrophage  phenotype  conversion,  and 
signaling)  as  the  main  driver  of  effective  continuation  and  reso¬ 
lution  of  inflammation  (61,  70,  71).  The  migrating  blood  mono¬ 
cytes  are  polarized  into  a  wide  spectrum  of  macrophage 
phenotypes,  and  the  two  ends  of  that  spectrum  can  be  character¬ 
ized  as  pro-inflammatory  macrophage  phenotypes  (similar  to  the 
“classically  activated,”  or  Ml,  phenotype  induced  in  vitro  by 
bacterial  LPS  and  IFN-y)  and  anti-inflammatory  phenotypes 
(similar  to  the  “alternatively  activated,”  or  M2,  phenotype  induced 
in  vitro  by  IL-4  and  IL-13),  respectively.  According  to  the  scenario 
modeled  in  this  work,  Ml -type  macrophages  differentiate  into 
M2-type  macrophages,  which  are  responsible  for  inflammation 
resolution  (59,  60,  70,  71).  However,  other  scenarios  of  macro¬ 
phage  phenotype  progression  are  conceivable  and  are  supported 
by  existing  data.  For  example,  in  mouse  wounds,  distinct  circu¬ 
lating  monocyte  subsets  can  be  recruited  to  the  injury  site  in  a 
temporally  ordered  manner  and  subsequently  differentiate  into 
distinct  macrophage  phenotypes  (59).  Specifically,  the  first  and 
second  groups  of  monocytes  sequentially  entering  the  wound  site 
give  rise  to  macrophages  with  phenotypes  resembling  those  for 
classical  and  alternative  activation,  respectively.  Interestingly,  this 
scenario  is  IL-4  and  IL-13  independent,  as  evidenced  by  the 
presence  of  M2-type  macrophages  in  the  wounds  of  IL-4Ra- 
knockout  mice  and  by  the  negligible  levels  of  IL-4  and  IL-13  at 
the  wound  site  (59,  72).  Furthermore,  in  some  cases,  accumulation 
of  substantial  numbers  of  M2-type  macrophages  is  driven  pri¬ 
marily  by  IL-4-dependent  proliferation  of  tissue-resident  macro¬ 
phages,  independently  of  blood  monocyte  recruitment  (73).  Such 
scenarios  potentially  can  be  stimulated  and  modulated  in  vivo  by 
a  myriad  of  molecular  and  cellular  processes.  Indeed,  it  was 
suggested  that  the  Ml -type  to  M2-type  phenotypic  transitions  are 
impacted  by  various  factors,  including  macrophage  interaction 
with  the  extracellular  matrix,  neutrophils,  cellular  debris,  and 
other  soluble  molecular  mediators,  such  as  IL-10  and  TGF-|3  (59, 
60,  74,  75). 

In  addition  to  local  regulation,  systemic  influences,  including 
immune  complexes,  PGs,  and  glucocorticoids,  may  determine  the 
phenotypic  fate  of  the  monocytes  migrating  to  the  inflammation 
site.  For  example,  macrophages  differentiated  in  the  presence  of 
the  glucocorticoid  dexamethasone  exhibit  increased  sensitivity  to 
TGF-J3  resulting  from  an  increase  in  TGF-J3  receptor  expression 
and  induced  lipid  uptake  (76).  However,  the  relative  contributions 
of  these  numerous  factors  to  macrophage  polarization,  as  well  as 
the  complex  and  synergistic  interactions  between  the  factors,  are 
unknown. 

The  choice  of  driver  for  the  macrophage  phenotypic  transition 
for  our  model  (i.e.,  pro-inflammatory  macrophage  phagocytosis  of 
apoptotic  neutrophils)  was  based  on  the  availability  of  experi¬ 
mental  data  that  could  be  translated  into  a  modeling  framework. 
Indeed,  this  mechanism,  which  is  supported  by  experimental  ev¬ 
idence  (40-42,  59,  75),  robustly  drove  the  macrophage  phenotypic 
transitions  in  our  simulations.  Although  IL-4  and  IL-13  might 
further  modulate  this  mechanism,  their  primary  contribution 
in  vivo  appears  to  be  the  connection  between  innate  and  adaptive 
immune  systems,  because  their  sustained  production  is  predomi¬ 
nantly  due  to  adaptive  immune  system  cells  (74).  However,  be¬ 
cause  of  our  emphasis  on  the  self-modulating  function  of  innate 
immunity  and  the  limitations  in  our  understanding  of  the  inter¬ 
actions  between  innate  and  adaptive  immunity,  we  do  not  model 
that  connection  in  this  study.  The  possible  contributions  of  other 
molecular  mediators  to  the  functional  effects  of  macrophage  po¬ 
larization  were  partially  reflected  through  the  IL-10-  and  TGF-(3- 


mediated  reduction  in  pro-inflammatory  cytokine  production  rates. 
Although  our  approach  to  modeling  macrophage  phenotype  con¬ 
version  is,  by  necessity,  an  oversimplification  of  the  molecular 
and  cellular  events  shaping  the  transition  to  the  inflammation- 
resolution  phase,  it  may  provide  a  reference  point  for  future 
studies  specifically  focused  on  macrophage-polarization  mecha¬ 
nisms. 

The  main  limitations  of  this  study  arise  from  the  simplifications 
needed  to  develop  a  mathematical  description  of  the  incredibly 
complex  cellular  and  molecular  interactions  that  occur  during 
normal  and  pathological  inflammation.  First,  our  model  focuses  on 
general  aspects  of  inflammation  and  does  not  reflect  the  specifics  of 
certain  inflammation  scenarios,  such  as  burns  or  bacterial  infection. 
Moreover,  we  do  not  explicitly  define  tissue-specific  cells,  such  as 
resident  macrophages,  mast  cells,  and  connective  tissue  cells.  Our 
model  is  based  on  the  notion  that  key  participants  (neutrophils  and 
macrophages)  and  interactions  in  the  local  inflammatory  process 
may  be  sufficiently  similar  across  many,  or  most,  inflammatory 
situations  (3,  11,  14,  16,  52,  57,  77,  78).  Therefore,  although  we 
modeled  chronic  inflammation  in  (traumatic)  wounds,  we  ex¬ 
pect  that  our  findings  may  be  applicable  in  other  contexts,  which 
is  supported  by  comparisons  of  our  modeling  results  with  exper¬ 
imental  studies  of  noninjury-induced  inflammation  (Figs.  3,  5). 
Second,  our  main  goal  was  to  develop  a  computational  model 
relevant  for  human  inflammation,  yet  our  model  training  and 
validation  were  performed  using  data  for  rat  and  mouse  (in  ad¬ 
dition  to  human)  experimental  models.  Although  interspecies 
differences  in  the  details  of  the  inflammatory  response  may  occur, 
the  use  of  animal  models  is  a  mainstay  of  immunity  research  (79). 
As  a  result,  several  experimental  datasets  that  were  indispensable 
for  our  modeling  efforts  had  been  obtained  with  animal  models, 
whereas  analogous  datasets  for  human  experimental  models  were 
not  available.  Third,  in  most  cases,  we  assume  constant  (but  cell 
phenotype  dependent)  cytokine  production  and  cell  activation 
rates  and  do  not  explicitly  consider  intracellular  signaling  path¬ 
ways  that  regulate  such  production  and  activation.  We  chose  this 
approach  because  separation  of  intra-  and  intercellular  kinetic 
scales  could  sharpen  our  focus  on  cell  phenotype  dynamics  and 
general  molecular  mediator-production  patterns,  which  aligned 
with  the  main  goals  of  this  work.  Fourth,  in  the  model,  we  reduced 
the  variety  of  macrophage  phenotypes  that  are  known  to  exist  at 
inflammatory  sites  to  only  two,  which  gave  us  the  benefit  of  mod¬ 
eling  tractability.  Finally,  our  model  represents  pro-inflammatory 
macrophage  phenotype  conversion  as  the  only  source  of  anti¬ 
inflammatory  macrophages,  whereas  independently  activated, 
“second- wave”  migration  of  anti-inflammatory  macrophages  into 
the  inflammation  site  is  a  plausible  alternative  (60,  71).  Further 
studies  are  needed  to  develop  and  validate  more  detailed  repre¬ 
sentations  of  the  inflammatory  process  and  its  perturbations. 

Because  of  the  recognition  of  the  role  of  chronic  inflammation 
in  different  pathologies,  the  need  for  efficient  therapies  and  early 
diagnostic  molecular  tools  for  chronic  inflammation  has  rapidly 
increased.  Our  modeling  suggests  that  narrowly  targeting  one 
single  mechanism  may  not  always  be  an  efficacious  strategy  to  treat 
chronic  inflammation,  because  it  can  be  triggered  and  maintained 
by  different  mechanistic  factors  at  different  time  points  (Fig.  4). 
This  might  explain  the  limited  efficacy  of  some  currently  used 
therapeutic  interventions,  such  as  nonsteroidal  anti-inflammatory 
drugs,  which  target  the  production  of  PGs  (80).  Likewise,  targeted 
neutralization  of  a  single  pro-inflammatory  cytokine,  such  as  IL- 
1  (3  or  TNF-a,  used  to  treat  chronic  inflammatory  pathologies  as¬ 
sociated  with  cancer  (81),  rheumatoid  arthritis  (82),  infection  (83), 
and  impaired  cutaneous  wounds  (84,  85),  has  seen  limited  success. 
Rigorous  experimental  testing  and  extensive  clinical  evaluation 
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of  our  computational  modeling  results  might  contribute  to  the  de¬ 
velopment  of  timely  therapies  for  chronic  inflammation  and  provide 
a  panel  of  informative  molecular  indicators  for  objective  clinical 
assessment  of  inflammatory  scenarios  in  injury  or  disease. 

Software  availability 

The  MATLAB  code  for  our  computational  analyses  is  freely 
available  and  can  be  downloaded  from  the  Journal  of  Immu¬ 
nology  Web  site. 
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